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Abstract 

Based on Density Functional Theory (DFT) with non-collinear-magnetism formu- 
lations, we have calculated the magnetism of single atom of 3d transition metals 
and the magnetic anisotropy of supported Ni chains on Au(110)-(lx2) surface. Our 
results show that surface relaxations enhance the orbital moments of left-end ele- 
ments (Ti,V) and quench the orbital moments of right-end elements (Fe,Co,Ni) on 
the Au(110)-(lx2) surface. This is because Ti and V atoms raise their positions 
above the trough of the reconstructed Au(llO) surface and Fe, Co, Ni atoms trap 
more deeply in the trough after the surface relaxations. From the study of magnetism 
of Ni one-dimensional atomic chains, their magnetic anisotropy are closely relate to 
their orbital quenching. The easy magnetized axis changes from the direction par- 
allel to the chains to the direction perpendicular to the Ni chains when they absorb 
on the surface. These one-dimensional Ni chains in the trough of the reconstructed 
Au(llO) surface have ferromagnetic order and are favored with chain-length less 
than 6. 

Key words: Non-Collinear Density Functional Theory, Surface Magnetism, 
one-dimensional Magnetism, Magnetic Crystalline Anisotropy 
PACS: 75.70.Ak, 75.30.Gw, 75.75.+a, 71.15.Mb 



1 Introduction 



The high-density magnetic recording and the memory devices require the 
strong magnetic anisotropy. The low-dimensional materials due to the reduced 
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dimensionality generally have a favorable direction and the physical prop- 
erties along this direction are generally different from the other directions. 
Thus we have the most possibilities to find the strong magnetic anisotropy 
materials from these low- dimensional materials. On the other hand, the re- 
duced dimensionality reduces the atomic coordination number and enhances 
the spin and orbital magnetic moments of the materials. Since the pioneer- 
ing experiments for the magnetism of Fe strips on the W(110) and Cu(lll) 
surfaces [1-4], the researches along this direction have been extended to oth- 
ers quasi-one-dimensional systems such as the monatomic magnetic chains on 
vicinal surface [5-7] and Co magnetic dot-chain on Ru(0001) surface [8]. The 
transition from in-plane to out-of-plane magnetic anisotropy had been found 
in Fe nano-structure when it is approach to one-dimensional limit [9]. Be- 
sides the celebrating properties such as the strong magnetic anisotropy, the 
enhanced spin and orbital moments, some new phenomena have been found, 
such as the temperature and time dependent magnetization [1,2] which make 
the realistic applications of these novel experimental devices face new serious 
problem such as the instability of magnetic structure. 

It is valuable to investigate theoretically the magnetic stability of surface- 
supported novel structures [10]. Many theoretical and computational methods 
have been used to studied the magnetism of one-dimensional structure on 
surface such as the KKR [11-14], and TB-LMTO [15] and PAW methods [16, 
17] based on Density Functional Theory, and the numerically tight-binding 
self-consistent calculation [18] and analytically tight-binding calculation [19]. 
Some common arguments have been found such as the enhancement of surface 
magnetism, the magnetic anisotropy are close related to the orbital moments 
and the spin-orbit coupling interaction (SOC). 

The orbital polarization [20-22] can remedy the small orbital moments in 
the calculations of density functional theory to approach the experimental 
values, and especially reproduce the correct phase transition in Ce companying 
with the volume collapses [21]. The orbital polarization plays the opposite 
role to crystal field which quenches the orbital moment to small value in 
crystal due to the splitting of ground state energy levels. However there are 
still many contradictions among the already known theoretical results. The 
density functional calculations with OP (orbital polarization) + SOC (spin- 
orbit coupling) have given the correct orbital moments but too large MAE 
(Magnetic Anisotropy Energy) for tetragonal and cubic Ni. If the OP term 
switches off the theoretical value of the magnetic anisotropy energy is close its 
experimental values [23]. The magnetic anisotropy energy of CoPt calculated 
using the so-called c-RPA [24] are weakly dependent on the new introduced OP 
interaction. The magnetic anisotropy energies are surprisingly consistent with 
experiments even only under LSDA approximation without the OP effects. 
The orbital moments of transition metals are well consistent with experiments. 
This method has been combined with PAW method to calculate quasi-particle 
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spectrum of materials [25]. 



Based on Brook's theory, the magnetic anisotropy is induced by the spin- 
orbital interaction [26]. Although the orbital moments are underestimated, 
it's still valuable to study the magnetic anisotropy properties in the absence 
of (OP) orbital polarization term. Theoretical works based on Density Func- 
tional Theory have also discovered others interesting properties of surface 
supported magnetic chains. Such as in Ref. [17], the authors have studied the 
magnetism of the ultrathin wires of fifth and sixth row elements supported 
on Cu(117) and Ag(117) vicinal surfaces. They found that only metals with 
a half-filled d band are found to have magnetic order on Cu(117) surface, ad- 
ditionally the ferromagnetic order is energetically stable. On the contrary, on 
the Ag(117) surface, the anti-ferromagnetic order is stable. The magnetism of 
mono-atomic chain is strongly dependent on their local environment absorbed 
on the surfaces. 

Besides the high-index vicinal surfaces such as Ag(117), the reconstructed 
Au(110)-(lx2) surface decorated with one-dimensional troughs along the closed- 
packed A-[110] direction can be used as a template to grow one-dimensional 
nano-structure. Experimentally, the distribution of lengths of Ni chains shows 
that most Ni chains are short than 6 [27]. Our molecular simulation addi- 
tionally shows that the lattice misfit between Ni chains and the Au substrate 
can induce the instability of Ni chains and makes them break into short seg- 
ments [28]. 

In this work we study the magnetism of mono-atomic Ni chains supported 
on Au(110)-(lx2) surface using the same theoretical method as in Ref. [17]. 
Our studies focus on the magnetic anisotropy of surface-supported magnetic 
atomic chains. We perform the non-collinear magnetic calculations including 
the spin-orbit coupling interaction but without orbital polarization interac- 
tion. We have calculated the magnetic anisotropy energies of the supported 
Ni chains and found an off-plane easy magnetization axis perpendicular to the 
chains and the surface. On the contrary, the easy magnetization axes of the 
free-standing Ni chains are parallel to these chains. The change of the easy 
magnetization axes of the surface-supported chains are also found in the oth- 
ers systems [18]. Our results also indicate the closed relationship between the 
magnetic anisotropy and the orbital magnetism. 

The paper is arranged as follow. In the next section, we introduce the theo- 
retical methods used in this paper. The third section will present the results 
for single 3d transition metal atoms absorbed on Au(110)-(lx2) and Au(lll) 
surface. The magnetism of a single magnetic atom on a surface is important 
to understand the magnetism of materials with more complex structures such 
as the one- dimensional nano-structure. The reduced dimensionality and coor- 
dination number for surface still enhance the magnetism of single absorbed 
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atoms [29-33]. The results of short one- dimensional mono-atomic Ni chains are 
included in the forth section. Our studies focus on the magnetic anisotropy of 
Ni mono-atomic chains. Finally, we conclude our results. 



2 Theoretical Methods: Non-collinear Density Functional Theory 
and PAW methods 

We have calculated magnetism of Ni n (n=l-5) chains supported on the Au(llO)- 
(1x2) surface based on Density Functional Theory [34,35] and the Methods 
of Projection of Augmentation Wave (PAW) [36] with the plane-wave base set 
and GGA Perdew-Burke-Ernzerhof's exchange-correlation potential [37]. The 
PAW methods used in VASP program [38-40] is as accurate as frozen-core all- 
electron methods and hopefully improved to include all-electron relaxation. 
Our calculations include the non-collinear effects and the spin-orbit coupling 
which is proved important to heavy metals such as gold. In non-collinear for- 
mula, a magnetic moment as a vector can point to any direction in space. The 
orbital polarization (OP) interaction isn't included in our calculations. 

The basic theory used in this paper is Hobbs's fully unconstrained non-collinear 
DFT formulas based on plane wave basis sets and all-electron PAW meth- 
ods [40,41]. Instead of distributing local quantization axes for every atoms, 
which is used in programs based on atom sphere approximation (ASA) or 
others analogous methods [42], Hobbs's theory adopts only one global quanti- 
zation axis and the vectors of magnetic electronic density are varying smoothly 
and continuously, which is more suited for the calculations of itinerant mag- 
netism in transition metal materials. 

In non-collinear density functional theory [40-42], the total energy is the func- 
tional of a density matrix n a p where a and f3 are the spin index along a defined 
quantization axis, a =| or j represent the spin angular momentum point to 
the positive direction or negative direction of the quantization axis. The sum 
of the diagonal elements is the total charge density, that is n Tr = J2a n aa- In 
a non-collinear magnetic system, the off-diagonal elements are non-zero. The 
2x2 matrix can be expanded using the complex matrix basis (/, cr x , a y) a z ) 
and expressed as 

n a( 3 = [n Tr S a( s + rh ■ a a/3 ]/2, (1) 

where m = J2 a /3 n Pa • is the magnetic electronic density and a is the Pauli 
matrix. In PAW methods, the pseudo-density matrix can be expressed using 
the pseudo-wave functions as 
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n 



(2) 



where /„ is the occupation number. The total electronic density matrix n a p = 
n a /3 + n ai3~ ^afii where nig and nig are the on-site electronic density matrixes. 
The total energy as the sum of three parts E = E + E 1 — E 1 , where E 
the smooth functional of the pseudo-density matrix n a p(r) is evaluated on a 
regular grid, E 1 and E l the functional of on-site electronic density matrixes 
n^p and n l a p are calculated on radical support grid around every atom. 

The Kohn-Sham Hamilton and Kohn-Sham equation can be obtained from 
the variation of total energy to the soft pseudo-density matrix and written as 



H°0[n] = -±AS°* + v$ f (3) 

(ij) 

J2H^\^ n )=e n S aa \K)- 

where S aa is the overlapping operator, Dff, 1 Dff and 1 D°f are the non- 
local interaction obtained by the PAW transformation, which are equivalent 
to the nonlocal pseudo-potential in the ultrasoft-potential methods. In non- 
collinear DFT formulas, the one-electron effective potential v^ff = VH(r)5 a p + 

Vxc{ r )8 a /3 + b(r) ■ <? a/3 , which is the summation of electrostatic potential vh(t), 
the nonmagnetic part of GGA exchange-correlation correction v xc (r) and the 

— * — * 

magnetic exchange-correlation potential b(r) ■ a a/j , where b = 5E xc /8fh(r). 
The magnetic exchange-correlation potential will potentially enhance mag- 
netic moments in non-collinear calculations. 

The spin-orbit coupling enters into PAW formulas by all-electron part of PAW 
Hamiltonian and is expressed as 



K — 1 dV -> 

where pi is the projector function and 0, all-electron wave function. The one- 
electron effective potential is modified v"^ — > v^ff + HgQ T , and the same 
time, one part of nonlocal PAW potentials is modified according to D°f — > 
Dff = J2l f(Veff( r ) + Hsoi( r ))Qij(. r )d r - The roles of spin-orbit interaction in 
and 1 D'*f are similar to , however they are calculated only within 
atom sphere. In the calculations of real material, they adopt the values in the 
calculations of isolate atom and don't update in Kohn-Sham self-consistent 
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loop. More information on non-collinear density functional theory within PAW 
methods can be found in the reference [40]. 

The atoms can have their initial magnetic moments by constructing the mag- 
netic electronic density and the magnetic moments around each of atoms are 
equal to their desired values. If the initial total magnetic moments point n 
direction, generally the direction of total moment changes in Kohn-Sham self- 
consistent loop. We can use the constrained DFT calculations to fix the di- 
rection of magnetic moments. Our self-consistent calculations show that the 
directions of magnetic moments change very small in ferromagnetic materials 
although using the unconstrained DFT calculations (see the table. (2) in sec- 
tion 4.1 ). In real ferromagnetic materials, the directions of magnetic moments 
are hard to change, once the magnetic domains have formed. This is because 
the change of magnetization direction must include the changes of magnetic 
moments of all magnetic atoms in the materials. An external magnetic field 
can force to change the direction of magnetic moments. Of course there exists 
an easy magnetization axis which is stable direction of magnetic moments. 
We need the 'strongest' magnetic field to change direction of magnetic mo- 
ments away from the direction of the easy magnetization axis. However for 
the materials with more complex magnetic structure such as the canted and 
spiral magnetic structures, the relaxations of directions of magnetic moments 
are intrinsic, the use of unconstrained DFT is essentially important [40,41]. 



3 The magnetism of 3d transition metal atoms absorbed on Au(llO)- 
(1x2) and Au(lll) surface 

The missing row Au(110) surface reconstructs from Au(110) surface by hav- 
ing missed closed-packed Au chains every others, and forms one-dimensional 
trough along the close-packed A- [110] direction. We construct a crystal slab 
with 4 atomic layers including 28 gold atoms. At the first step we optimized 
the lattice constant of gold crystal which is about 4.2A. In the procession of 
optimization, the RMM-DIIS algorithm [38] is applied to the total-energy min- 
imum with plane-wave energy cutoff 229. 9eV and the 6x6x6 Monkhorst-Pack 
K points mesh. Now the size of super-cell is 11.8Ax8.4Ax5.94A 

We add a vacuum layer by extending the super-cell along C-[110] direction. 
The size along C-[110] direction of super-cell including vacuum layer is equal to 
16A The lattice constant of surface is generally smaller than that of the crys- 
tal. Thus we optimized the slab with vacuum layer together. Based on the op- 
timized lattice constant, the size of the super-cell is 11.523Ax8.148Ax 15.52A 
The vacuum layer is about twice thicker than thickness of the slab. For the 
optimized Au(lll) surface, the size of super-cell is about 11.5Ax9.97Ax20A 
including 64 gold atoms. In above optimizations, we only change the lattice 
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constant. The initial structures are constructed by placing the 3d transition 
metal atoms on the hollow sites of the Au(lll) and Au(110)-(lx2) surface 
in the trough (Fig. 1). The initial structures are optimized for all atoms ex- 
cept for the bottom atoms by having used the conjugate-gradient Methods. 
The K-points meshes of 2x2x1 or T point are used to sampling the first 
Brillouin-zone in the corresponding calculations of the electronic structure for 
Au(110)-(lx2) and Au(lll) surfaces respectively. In this work we use smaller 
k-points mesh to more efficiently optimize surface structures. 

We use both collinear and non-collinear density functional calculations but 
mainly present the non-collinear results. In non-collinear DFT calculation, we 
use the optimized structures from collinear DFT calculation. Based on the 
relaxed structures, we calculate more exactly the electronic structures using 
the RMM-DIIS algorithm [38] with the (6x6x1) Monkhorst-Pack grids sam- 
pling the first Brillouin zone. If the changes of the total energies are smaller 
than O.OOOleV between two electronic self-consistent (SC) steps the SC-loops 
break. We set the Methfessel-Paxton smearing width equal to 0.20eV to accel- 
erate the speed of convergence. The energy-cutoffs of plane-wave for different 
atoms are summarized in table (1). The magnetic moments for a single atom 
are calculated by considering the electrons in a Wigner-Seitz sphere centered 
at the position of the atom. The Wigner Seitz Radii for 3d transition metal 
atoms are also summarized in table (1). We choose the C-[110] direction as 
the quantization axis. We emphasize the influence of surface deformation on 
the orbital moments of the absorbed atoms. The orbital polarization interac- 
tion omits in this work but includes the spin-orbit coupling interaction. We 
obtain smaller orbital moments compared with pervious density functional 
calculation including the orbital polarization interaction [44]. 

3.1 The influence of structural relaxations on the magnetism of a single 
absorbed 3d transition atom 

It's seem that collinear DFT calculation is enough to the study of magnetism 
of a single transition metal absorbed on the gold surface. However, if we want 
to know the direction of magnetic moments relative to the surface, we must use 
the non-collinear DFT formulas. These results for single absorbed atom will 
be compared with the results of one-dimensional atomic chains. At first step, 
we calculate the magnetic moments of the absorbed 3d atoms on the Au(llO)- 
(1x2) and Au(lll) surfaces including the spin-orbit coupling interactions and 
non-collinear calculations. From Fig. 2, we can see that the spin moments 
reach the maximum in the middle of the group. The large spin moment of the 
absorbed Mn atom is about 4.163/i s on the Au(110)-(lx2) surface. We can 
see from the left panel of the figure that the spin moment of the absorbed Cr 
atom is the second largest when including the spin-orbit coupling interactions 
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and non-collinear calculations. The change of magnetic moments per absorbed 
magnetic atom across the 3d row has been found in many other ab initio 
calculations such as 3d atoms on Au(OOl), Ag(OOl) [45] and Cu(OOl) [46] 
surfaces. Our results are consistent very well with these results although the 
energy cutoffs of plane-waves are not too large in our calculations. 

Both the electronic correlations and the crystal field have significant influence 
on the orbital moments of the absorbed 3d atoms. If the electronic correlation 
is stronger than the crystal fields, the orbital moment is large, otherwise small. 
The absorbed Ti,V,Co atoms have visible orbital moments 0.086/z B , 0.082/i B 
and 0.092// B on the Au(110)-(lx2) surface, and0.006/i s , 0.071// B and 0.247// B 
on the Au(lll) surface respectively. Cabria, et.al obtained larger orbital mo- 
ment (about 0.5// s ) of Fe and Co atoms absorbed on Au(OOl) surface using 
the spin polarization relativistic KKR methods including the orbital polariza- 
tion term [44,45]. The usual exchange-correlation potentials (LDA or GGA) 
underestimate the electronic correlations such as the Coulomb correlation and 
the orbital polarization, thus the obtained orbital moments in this work are 
generally small compared with experimental values [20-22, 44] although our 
work includes the spin-orbit coupling interaction. The small values of orbital 
moments of the absorbed Cr atoms in our calculations are not related to the 
crystal field but to its electronic structure. In the individual Cr atom, five 
3d-electrons half-fill the 3d states, the total orbital moment is very small. 

The coordination number of surface atoms is generally smaller than those in 
bulk. So the magnetism of surface is generally stronger than that in bulk. The 
supported atoms on the surface are in the environment similar to the surface 
atoms. Based on the same logic, the supported atoms on the surface have pos- 
sibly strong magnetism. The lack of the orbital polarization isn't the obstacle 
to study the structural influence on the orbital moments. In order to clearly 
illustrate the effects of the crystal deformations (or the cubic distortions), we 
compare the results on the deformed surface with that on the prefect Au(110)- 
(1x2) surface. The absorbed 3d atoms on the prefect surface still modify their 
positions to reach their stable positions, while the surface atoms are fixed. 
We find from Fig. 3 that the surface relaxations generally quench the orbital 
moments of the absorbed 3d atoms except for the absorbed V and Ti atoms. 
The changes of the orbital moments in response to the relaxations are closely 
related to the changes of the depth of the 3d atoms embedded in the trough 
of the reconstructed Au(110) surface. The absorbed atoms are deeper in the 
through, their orbital moments are quenched to smaller values due to the 
stronger crystal field. The absorbed V atom rises above the top row of the 
trough after the relaxations and its orbital moment enhances, opposes to the 
orbital quenching for the other 3d atoms except for the absorbed Ti atom. This 
is due to the weaker crystal field above the surface than that in the trough. 
The absorbed Ti atom has almost the same height as the top row after the 
relaxations and its orbital moment slightly increases. Thus our results indicate 
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that the surface relaxation decreases the orbital moments of the absorbed 3d 
atoms with the excess half-filled 3d states and increases the orbital moments 
of the absorbed atoms with less half-filled 3d states on Au(110)-(lx2) surface. 

The 3d density of states of the absorbed 3d atoms on the two surfaces are 
shown in Fig. 5. The absorbed Ti, V, Co atoms on the Au(110)-(lx2) surface 
and V, Co atoms on the Au(lll) surface have large 3d densities of states at the 
Fermi energies and large orbital moments too. We get the same arguments on 
the magnetism of Ni n (n=l-5) chains. Large density of state N(Ef) at Fermi 
energy is advantaged to the formation of the orbital moments. The single Ni 
atom has small magnetic moment, however we can see in section 4 the Ni 
clusters such as one-dimensional Ni chains have significant magnetic moments 
on the Au surface. The exchange interactions between Ni atoms prevent from 
the decay of magnetic moments of Ni atoms. 



4 The magnetism of Ni atomic chains 

The initial Ni chains are located in the trough along the A- [110] direction with 
nearest-nearby distance about 2.88A the same as the nearest-nearby distance 
of the missed Au chain [Fig. 1]. The surface slab includes 4 atomic layers 
with 56 Au atoms and the size of the super-cell is 23.046Ax8.148Ax 15.52A. 
We have relaxed the initial structures using the conjugate-gradient Methods. 
The 2x2x1 K-points mesh is used to sampling Brillouin-zone in the corre- 
sponding electronic-structure calculations. Based on the relaxed structures, 
we calculated their electronic structures more exactly using RMM-DIIS algo- 
rithm with the (4x4x1) Monkhorst-Pack grids sampling the first Brillouin 
zone. The other information is the same as in the calculations for single 3d 
atoms. 

We have also calculated the magnetism of free-standing Ni chains by removing 
all surface atoms but still preserving in the same box cell. There are mirror 
atomic chains because of periodic boundary condition along three directions. 
The interaction of free-standing chains with its mirror atomic chains are small 
and negligible if the box cell is large enough. We have only chosen short chains 
Ni n where n less than 6 and large box size 23A along the direction parallel the 
chain so that the ends of chains are far away from the boundaries of the box 
cell. The energetically favorable configuration of one-dimensional free-standing 
chains has zigzag shapes and not straight line. The Ni n atomic chains in the 
trough of Au(110)-(lx2) surface are straight. As we optimize the structures of 
free-standing chains, all atoms are only allowed modified their position along 
the direction of chains and the optimized free-standing Ni chains still keep a 
straight line shape. We try to compare the results of surface-support straight 
chains with those of free-standing straight chains. The quantization axis is 
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still along the C-[110] direction in this section. 



4-1 Magnetism of Ni n (n—1-5) 

All Ni atoms in the cell have their initial magnetic moments in our calcula- 
tions. We calculate the total energies, the spin and orbital magnetic moments 
of the supported Ni chains when they are magnetized along A- [110] direction 
parallel to the chains or C-[110] direction perpendicular to the chains respec- 
tively. The magnetization of an atomic chain in this paper means that the 
magnetic moments for every Ni atom align along the same direction. The val- 
ues of initial moments of Ni atoms are all the same, however the values after 
the self-consistent calculations are probably not the same for all Ni atoms. The 
orientations of magnetic moment may manually change to simulate the dif- 
ferent magnetizing directions. In a single self-consistent calculation, the direc- 
tions of magnetic moments change very small for the ferromagnetic Ni chains 
although we have used the unconstrained Density Functional calculations (see 
the table. (2)). 

The chains become short after the relaxations and the average nearest-nearby 
distance of Ni 5 chain is about 2.75 A shorter than 2.88 A that of initial struc- 
ture [the insert figure of Fig. 5(b) and Fig.l]. We choose two magnetizing direc- 
tions: A- [110] direction parallel to the chain and C-[110] direction perpendic- 
ular to the chain. The size-dependence magnetic moments are shown in Fig. 5 
(a). From this figure we can see that the spin moments of a single absorbed 
Ni atom are small, 0.147 /a B for the magnetization along C-[110] direction and 
0.127// s for A-[ll0] direction. The spin moments pre atom of the Ni 2 chain 
increase to 0.455 li b for A-[110] direction and 0.464 \x B for C-[110] directions. 
The spin moments of Ni 3 have slightly decrease for C-[110] direction but large 
decrease for A- [110] direction. The spin moments decrease slightly and try 
to keep a constant as the chain length increases above 4. Our results show 
that the spin moments of Ni atoms on the reconstructed Au(110) surface are 
smaller than the value 0.675/zs in crystal Ni having been calculated using the 
same methods in our work. 

The magnetic order is an important aspect of material magnetism. Experimen- 
tally the long ferromagnetic order is found for Co chains stabilized above a 
finite temperature on Pt(997) surface [5]. Above the temperature, the long fer- 
romagnetic order is destroyed and changed into the short ferromagnetic frag- 
ments. In our calculations, the short Ni chains are ferromagnetic. The Fig. 6 
shows the ferromagnetic order of the surface supported M5 chains for the mag- 
netization for C-[110] and A- [110] directions. The spin-polarizations of surface 
Gold atoms near the Ni chain are induced by the magnetic Ni atoms, which 
almost disappear for Gold atoms far from the Ni chain. We have also calcu- 
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lated the anti-ferromagnetic configurations, non-magnetic configurations and 
other non-collinear magnetic configurations and found their energies higher 
than that of the ferromagnetic configuration. Thus our calculations show that 
the Ni chains with ferromagnetic order are most energetically favorable. It is 
interesting to find from this figure that, when magnetizing the Ni chains par- 
allel to the chain, the orbital moments are larger than that perpendicular to 
the chains. The spin moments for every Ni atom of the ferromagnetic Ni chain 
almost don't change their directions after the Kohn-Sham self-consistent cal- 
culations for all three magnetization directions [table. (2)]. The spin moment 
of the middle atom is smaller than that of other atoms and thus the short Ni 
chain isn't in prefect ferromagnetic order. The larger spin-moments for atoms 
at two ends are because of the smaller coordination number compared with 
the middle atoms. 



The interaction between supported Ni n chain and surface can be measured 
with the interaction energy Ej = E tot — E s — E' where E tot is the total 
energy, E s the energy of slab with relaxed surface and E* the energy of free- 
standing chain. We generally calculate the interaction energy per atom, which 
removes the effects of simple summation. The interaction energies are gener- 
ally negative. Fig. 5(b) shows that the interaction energies of the Ni chains 
with the surface increase with the chain lengths. In our classical molecular 
dynamics simulation, we have found the oscillation of interaction energy for 
longer chains induced by lattice misfit [28]. These results indicate that short 
Ni chains are more energetically favorable on Au(110)-(lx2) surface. This is 
consistent with the experiments [27] in which most of Ni chains are shorter 
than 6. Additionally, we calculate the density of states of the Ni n (n=l-5) 
chains. From the Fig. 7, we can see that the densities of states of 3d states 
of Ni 5 chain changes very small when the magnetization changes from C-[110] 
direction to A- [110] direction and the spin magnetic moment does also change 
very small [Fig. 7]. From the figure 5, the changes of orbital moments are 
more significant if the magnetizing directions change. Although the orbital 
moments are underestimated due to the absence of (OP) orbital polarization 
term, our results still show that magnetic anisotropy are closely related to 
the orbital degree of freedom. For compact magnetic islands on metal surface 
such as Co/Pt(lll) [47], the orbital moments decrease with increasing island 
sizes. This is because only atoms at the edge of island have significant contri- 
butions to the total orbital moments. However, for one-dimensional island all 
atoms have significant contribution to the total orbital moments. The orbital 
moments per atom keep constant as the lengths of chains beyond 4. 
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4-2 The magnetic anisotropy of N% chain 



In the above subsection we have chosen two magnetization directions A- [110] 
and C-[110]. In this subsection, we study the magnetic anisotropy of the 
NI5 chain and show that C-[110] direction is along the easy magnetization 
axis. The magnetic anisotropy energies in Si and S2 planes are defined as 
5E 1 = E c - E A and 5E 2 = E c - E B respectively. E A , E B and E c are the 
total energies when magnetizing the chain along A-[110] (9 = 0°), B-[001] 
(0 = 0°)and C-[110] (0 = 90°) directions [53]. The magnetic anisotropy of 
one-dimensional structure is more prominent than compact structures such as 
surface and bulk. Fig. 8 shows the total energies, the spin and orbital moments 
change with two magnetization angles in above two planes. We find from 
Fig. 8(a,b) that there are weak magnetic anisotropy in Si plane and relative 
strong the magnetic anisotropy in S 2 plane. The magnetic anisotropy energies 
are 0.158(meV/atom) and 0.524(meV/atom) respectively in the two planes. 
The shape anisotropy determines the large difference of magnetic anisotropy 
energy in two plane because of one-dimension characteristic of mono-atomic 
chain. In the same plane such as the plane Si, the magnetic anisotropy is 
determined by spin-orbit coupling. The easy magnetization directions in both 
Si and S 2 planes are all along C-[110] direction (9 = 90° and = 0°), which is 
perpendicular to the chain and the surface. The magnetic anisotropy energies 
can be represented as Kq + K\sin 2 {9) in Si plane and K' + K[sin 2 (<f)) in S 2 
plane [26,52]. The parameters A" =0.158 (meV), K\ =-0.158 (meV) in Si and 
A"q=0.0 (meV), A^=0.62 (meV) in S 2 are obtained by fitting the magnetic 
anisotropy curves using above two functions. The magnetic anisotropy isn't 
good fit with Kq + K\sin 2 {9) in Si. The underlaying reasons are (1) there 
isn't simple cubic symmetry in the region where the one-dimensional atomic 
chain absorbed; (2) the heterogeneously distributed the nearest nearby dis- 
tances along the chain induced by the lattice misfit between the chain and the 
substrate. On the other hand, the value of K x is close to the lowest limit of 
allowed error in total energy about ±0.1meV, thus more accurate calculations 
are required. Our results also show that A"i (or K[) isn't equal to the simple 
5E 1 = E c — E A (or 5E 2 = E c - E B ). 

Experimentally, the values of magnetic anisotropy energy are from 1 meV 
to 10 meV for magnetic clusters on metal surface and from 0.1 mV to 1.0 
meV for magnetic film. Based on our calculations, the magnetic anisotropy 
energies for Ni chains on the reconstructed Au(110) surface are close to the 
values of magnetic film. This is because the atoms on Ni chain replace the 
atoms of first layer of original surface and have almost the same coordinate 
number in their monolayer film. The experiment in Co/Pt(997) has shown 
that the anisotropy energy about 2 meV for the monatomic Co chains and 
0.14 meV for the monolayer Co film on Pt(997) surface [5]. From Fig. 8(e) 
we can find that the orbital moments reach the minimum when magnetizing 
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the Ni chain along the easy magnetization direction. Fig. 8(f) also shows that 
the orbital moments almost keep a constant when the magnetization angle 
changes in S2 plane. The directions of orbital moments for every Ni atoms on 
the chain are almost parallel to their spin moments when the magnetizations 
are along A- [110] and C-[110] axes, or anti-parallel to their spin moments for 
the magnetization along B-[001] axis. From the curve of magnetic anisotropy 
energy in S2 plane in Fig. 8(b) it's most energetically unfavorable when the 
spin moment of the Ni chain points to B-[001] direction (0 = 90°). 



4-3 The magnetic anisotropy of free-standing Ni 5 chain and Discussion 

In order to illustrate the surface influence on the magnetic anisotropy, as a 
comparison, we also calculate the magnetic anisotropy of free-standing Ni 5 
chain in the same box. All calculating details are the same as those of the 
surface-supported Nis chain except for removing the surface atoms. The av- 
erage nearest-neighbor distance is about 2.196 A shorter than 2.77 A for the 
surface-support Ni 5 chain and 2.5 A in Ni crystal. The changes of densities of 
states of 3d states for free-standing Ni chain are that the majority-spin-up den- 
sities of states are far from to the Fermi Energy. The larger densities of states 
at Fermi energy contributed by the minority-spin-down DOS indicate that 
there are larger spin and orbital moment magnetic moments for free-standing 
Ni chain. 

We only consider the changes of the total energies with the angle 0. From 
the Fig. 8 (g) we can see that the easy magnetization direction is now par- 
allel to the chain which is different from perpendicular to the chain for the 
surface-supported Ni 5 chain. For free-standing chain we find K =0.0 (meV) 
and Ki—SA (meV) in Si. The anisotropy energy is large compared with the 
surface-supported Ni 5 chain. The change of the sign of K\ compared with 
surface-supported chain means that the easy magnetization axis change in Si 
plane. The change of the easy magnetization axis had been found in other the- 
oretical works such as Co chains on Pd(110) surface [18]. Experimentally, the 
changes of easy magnetization axis are found in the growth of magnetic film on 
metal surface when the thickness of film is beyond a critical value d c . The tran- 
sitions from in-pane to perpendicular plane for Ni/Cu(100) [48] and inversely 
from perpendicular to in-plane for Fe/Ag(100) [49] and Co/Au(lll) [50] are 
found at critical coverage about 5-7 ML. The transitions are close related 
with the transitions from films with heterogeneously distributed strains to 
films with homogenously distributed strains once beyond critical thickness. 
The atomic chains in this work are found in initial stage of film growth and 
have the easy magnetization axis perpendicular to the surface. Whether the 
easy magnetization axis changes after a monolayer had completely grown on 
the surface is interesting topic in our following work. 
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The spin moments keep almost constant when the magnetization angle changes 
from 9 = 0° to 9 = 90°. The spin moment about 1.14//B per atom is signifi- 
cantly large than its values in Ni crystal 0.675/i^. The behavior of the orbital 
moment is similar to the surface-supported Ni 5 chain, that is, reaches the min- 
imum when perpendicularly magnetizing the chain. The underlaying reason 
is geometric, that is, the orbital currents (or Molecular Currents) of neighbor 
atoms have large overlapping and their cancelations are also large which is dis- 
advantageous to form larger orbital moments. The values of orbital moments 
from 0.2/ie to 0.6/zb are larger one power of the values of surface supported 
Ni atoms. The easy magnetization axis of free-standing Ni chain A- [110] is 
different from the easy magnetization axis of surface-supported chain C-[110]. 
The orbital moment reaches its maximum along the easy magnetization axis 
A-[110], 9 = 0° for free standing chain, which is consistent with the results for 
magnetic monolayer, [22,51] but reaches the minimum along the easy magne- 
tization C-[110] axis 9 = 90° for the surface-supported chain. Thus our results 
show that the surface effects change the easy magnetization axis and decrease 
both the magnetic anisotropy energy and the magnetic moment of free stand- 
ing chains because of the increase of coordination number for every atom on 
a surface. 



5 Conclusion 



Based on the Density Functional Theory, we have calculated the magnetic 
anisotropy of the supported Ni chains and free-standing chain. The easy mag- 
netization axis of the supported Ni chain is perpendicular to the chain and 
surface but parallel to the chain for free standing chain. The different easy 
magnetization axis for the supported and free standing chain is also found 
in density functional calculations of other system. Our results indicate the 
closed relationship between the magnetic anisotropy and orbital moments of 
the Ni chains. In the environments of crystal, the orbital moment is generally 
quenched to small value due to the splitting of crystal field. Compared with 
spin moment, orbital moment is more sensitive to the change of magnetiza- 
tion direction. The orbital quenching changes with the change of magnetized 
directions. The magnetic anisotropy of materials are determined both by the 
orbital quenching and spin-orbit coupling. For a single absorbed atom, the sur- 
face relaxations generally deform the surface and modify the positions of the 
absorbed atoms. If the relaxations lift the absorbed atom out of the trough on 
the constructed Au(110) surface, the orbital moments enhance due to weaker 
crystal fields and if the relaxations make the absorbed atoms deeper embed- 
ded in the trough the orbital moments quench to smaller values due to the 
stronger crystal fields. 
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Table 1 



The Wigner Seitz Radius- WSR (A) and Energy-Cutoff of plane wave-Ecut (eV) for 
PAW potentials of 3d transition metal T 3(i =V, Ti, Cr, Mn, Fe, Co, Ni, Cu. [43] The 
real energy cutoff for T3^/Au is Ecut(T3^/Au)= Max [Ecut(T3^),Ecut(Au)]. 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Au 


WSR 
Ecut 


1.323 
178.4 


1.323 
192.6 


1.323 
227.1 


1.323 
269.9 


1.302 
267.9 


1.302 
268.0 


1.286 
269.6 


1.312 
273.2 


1.503 
229.9 



Table 2 

The spin and orbital magnetic moments (Bohr) for every atoms of Nis chain when 
the magnetization directions are along A-[110] ,B-[001] and C-[110] respectively. 
This table shows that the directions of spin magnetic moments are almost unchanged 
in Kohn-Sham self-consistent loop although we have used the unconstrained DFT 
calculations. 





A- [110] 


B-[001] 


C-[110] 




M A 


M B 


M c 


M A 


M B 


M c 


Ma 


M B 


M c 


Nil 


0.434 


0.0 


-0.006 


0.0 


0.411 


0.0 


0.017 


0.0 


0.430 


Ni 2 


0.449 


0.0 


0.006 


0.0 


0.426 


0.0 


-0.012 


0.0 


0.449 


Ni 3 


0.352 


0.0 


0.000 


0.0 


0.323 


0.0 


0.001 


0.0 


0.366 


Ni 4 


0.449 


0.0 


-0.007 


0.0 


0.452 


0.0 


0.015 


0.0 


0.452 


Ni 5 


0.437 


0.0 


0.005 


0.0 


0.446 


0.0 


-0.015 


0.0 


0.432 




La 


L_B 


Lc 


La 


Lb 


Lc 


La 


L B 


Lc 


Nil 


0.105 


0.0 


0.000 


0.0 


-0.085 


0.0 


0.009 


0.0 


0.080 


Ni 2 


0.097 


0.0 


-0.010 


0.0 


-0.082 


0.0 


-0.012 


0.0 


0.093 


Ni 3 


0.092 


0.0 


0.000 


0.0 


-0.065 


0.0 


0.001 


0.0 


0.060 


Ni 4 


0.097 


0.0 


0.010 


0.0 


-0.088 


0.0 


-0.012 


0.0 


0.096 


Ni 5 


0.105 


0.0 


-0.001 


0.0 


-0.092 


0.0 


-0.009 


0.0 


0.079 
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Fig. 1. The schematic diagram of a Ni chain supported on the Au(110)-(lx2) surface. 
The trough of the reconstructed surface is along the A- [110] direction. The red 
spheres are the absorbed Ni atoms in the trough. The yellow spheres are the Au 
atoms of the surface. The plane including both C-[110] axis and A-[ll0] axis is 
named as Si, the plane including both C-[110] axis and B-[001] axis as S2. The 
rotating angle from any vector in Si to A-[110] axis labels as and the rotating 
angle from any vector in S2 to C-[110] axis labels as <p. 9 and <p are the magnetization 
angles in the two planes. 
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Fig. 3. The effects of structure-relaxation on orbital moments of the absorbed 3d 
transition metal atoms, (a). The heights of top-row and absorbed 3d atoms be- 
fore the full-structure relaxations. The stable positions are different for different 
3d atoms, which are obtained by the single- atom optimizations on the unrelaxed 
surface, (b). The heights of top row and absorbed 3d atoms after having done ful- 
l-structure relaxations. The V atom is higher above the top row after the full-relax- 
ation, (c). The schematic figure defines the heights of top row and absorbed atoms 
using two arrows pointing to the bottom of the trough, (d). The orbital moments 
change before and after the full-relaxations. 
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Fig. 4. The 3d densities of state (DOS) for the 3d atoms absorbed on Au(110)-(lx2) 
(solid lines) and Au(lll) (dot lines ) surface. 
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Fig. 5. (a) The spin ^M, and orbital A L, c h moments change with the length 
of Ni chains. and C M aren't the A and B components of the spin moment M. 
They are the spin magnetic moments when the magnetization along the A- [110] 
and C-[110] direction and are obtained from different self-consistent calculations. 
A L and A L have the same meaning for orbital moments, (b) Interaction energies 
between the Ni chains and the surface. The larger absolute values indicate the 
stronger interaction with the surface. 
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Fig. 6. The vector graph of the spin magnetic moments of the Ni chain which illus- 
trates the ferromagnetic order of the M5 chain based on our calculations. Especially, 
(a) and (b) are corresponding to the magnetization along C-[110] and A- [110] direc- 
tion respectively. The smaller spin magnetic moments of Au atoms and the orbital 
magnetic moments of Ni atoms are multiplied by 40 and 2 respectively to show 
them with the larger spin moments of Ni atoms together in the same figure clearly. 



24 




Energy (eV) 

Fig. 7. The projected DOS of 3d states of the supported N15 chain along the di- 
rections of magnetization when the chain is magnetized along different directions 
respectively (the solid lines for the C-[110] direction and the circles for the A-[110] 
direction). The spin- up (or spin-down) means along the positive (or negative) di- 
rection of two magnetization directions. This shows that spin degree of freedom is 
less relevant to magneto-crystalline anisotropy. The dots line is the DOS for free-s- 
tanding Ni chain when magnetizing it along C-[110] direction. 
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Fig. 8. The anisotropy energies (a,b), the spin (c,d) and orbital moment (e,f) of the 
supported N15 chain change with the magnetization angles 9 and 4> in the planes 
Si and S2 defined in above. The easy magnetization direction is perpendicular to 
the surface 6 = 90° and <f> = 0°. The anisotropy energies (g), spin moment (h) and 
orbital moment (i) of the free standing N15 chain change with the magnetization 
angles 9 in the planes Si. The easy magnetization direction is parallel to the chain 
6 = 0°. Error bars in (a,b) represent the lowest limit of allowed error in total energy 
about ztO.lmeV. The magnetic anisotropy energy curves are fitting using function 
Kq + K\sin 2 {6) for (a,g) and Kq + Kisin 2 ((j)) for (b)(red dash- lines). 
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